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LAPLACE TFjiflSgQKAiS 

Pirst Lecture 

Textbook: Gardner and Barnes, ’’Transients in Linear Systems” 

Wiley, 1942 

Other books by Berg, Bush, Cohen, Carson, Kurtz and Corcoran, 
Skilling, Stephens, and Churchill may be consulted. 

Aim of course: to provide a technique for solving practical 
transient problems. 

The classical method will be studied first, using as an ex- 
ample a series-peaked video amplifier stage. This problem will 
later be solved operationally. We assume a ’’unit step” current 

applied to the interstage coupling circuit, defined by: 

, t < O 

i -io , t > 0 

The circuit and terminology are as shown below: 




The problem will be solved when we have found e (t). The 

standard classical procedure would be to set uu a differential 

equation in which e (t) is the unknown and i (t) the driving func- 
e \/t n 

tion. However,. in this relatively simple problem, the driving func- 
tion is not differentiable, and the procedure described would fail. 
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Hence, even in this case we must use an artifice in order to arrive 
at e (t)* which consists of finding i2 (t) first, and then using the 
solution of the sub- problem to find e (t). The fundamental equations 



of the circuit may be written as follows: 
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In order to eliminate i /t we differentiate (2): 
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Rearranging this equation, we have for the differential equation 
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This may be put into a standard form by substituting 
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The differential equation is then 
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The solution of which is: 
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It is seen that the condition for oscillations i; 
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The constants and C2 may be evaluated from the initial conditions: 
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We have then: 

C, ■+ c^t i 0 - o (&) 

C, + C* («.- Vo.'-i'-) - o 

from which we have: 




The voltage e ( t } may then be hound from i 2 *' 




Combining terms, we have: 




This is the solution of the problem. There are a few general- 



izations that may be made about the classical method. It is found 
that in all networks in which the components are linear and inde- 
pendent of time, the differential equations are of the type known 
as linear equations with constant coefficients. Their solution 
consists of the complementary function, which is the general solution 
of the homogeneous differential equation obtained by putting the 
driving-force equal to zero, and a particular integral which depends 
on the form of the driving-force function. These are known to the 
engineer as the transient and steady-state responses respectively. 
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The transient part of the solution contains a sufficient number of 
arbitrary constants to enable it to meet all the prescribed initial 
conditions. The functions appearing in the transient are always 
exponential. If these exponentials are complex they always occur 
in complex conjugate pairs, so that their sum represents a real, 
exponentially damped oscillation. The difference between two 
real exponential functions occurs often and is called a surge function 
It has the form illustrated below: 

f(t) = e -e 





The two modes for a simple LRC circuit vary with increasing losses 
as shown below: 

S 




Pt 



With no losses, ( Q'Jalfeo ) the modes are at ±7 uj 0 , and re- 
present an undamped sine wave, since ^ + d J ^ co 0 X 

As losses increase, the modes move along the arc of a circle, re- 
presenting an exponentially decaying oscillation of lower frequency 
than UP? . Critical damping occurs when the modes coincide at 
££ - z $^7- - nt , and at this point there is only a single ex- 
ponential term ^-^ot # As the losses increase further, the modes 




separate, and move along the real axis, where at infinite loss one 
approaches ( “ ) and the other approaches zero. 

The distinction between transient and steady-state is very 
great in the classical method of solution, but it is a mathematical 
rather than physical distinction, which is lost in the operational 




In the Heaviside method, we obtain an operational equation from 

„ df , , nt 

the differential equation by substituting f- ^ <? ~ - U yjf 



This transforms an 
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integro-diff erential equation into an algeb^ic one. If this is 
solved for the unknown treating p as an algebraic quantity, the re- 
sult is a symbolic solution involving p and the driving function. 
This can be converted into a real time function if we know how to 
interpret the operational equations; we will now consider how this 
is done. The fundamental type of driving function is the unit step 
which is defined by: r „ +. y 

1 - f 0 5 

' r l l 9 t- >0 J 

Heaviside found, by working backwards from known solutions of simple 

differential equations, that the following interpretations of 

operational expressions could be made: 
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and, in general, we have: 
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It should be clearly understood what this type of equation 
means; it is not valid to cancel out the / sign, because the left 
hand side of the equation is not a product, but rather represents 
an operation on J . The right hand side, however, is an ordinary 
product, and simply means that the time function is cut off before 
t equals 0, which amounts to assuming zero for all initial con- 
ditions, We see that if some function F(?) is given, the result 
of operating with this function on the unit step "j may be found 
by expanding I (P) in a power series in -L . after which each term 
may be interpreted according to equation (ll). Heaviside called 
•this process "algebrizing" . By means of it, an answer in the form 
of a power series in t can be obtained for any F (P) that is analytic 
in ~~r . Some of these power series are recognized as well-known 



functions; in this way the following identities have been established 
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We see that a time function can also be obtained in the form 
of a series of exponential terms if we expand F (P) into partial 

u 

fractions of the form ^ and use the appropriate formula above. 
The general rule for this process is called the Heaviside expansion 
theorem. If F (P) is expressible as the ratio AJ(f'> /oif) of two 
rational algebraic polynomials, we have: 
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where the are the" h roots of the equation D (f) -d , and 

{jttvC) is the derivative evaluated at (T-T* 

it is seen that the different "modes" of network response appear 



evaluated at ff -(7^. 



as different solutions of D(f) * o , which ensures that all com- 
plex roots will appear in complex conjugate pairs. 

If the time function operated on by F (P) is not the unit 
step, it is usually easiest to express that time function in terms 
of a function G (P) operating on the unit step, so that the pro- 
duct F (P) G { P ) then becomes the function to interpret. For 



example , 
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Second Lecture 



The last lecture covered classical and Heaviside methods of 
solving differential equations. 

The reason for success of Heaviside method was that it re- 
duced the differential equation to an algebraic one which could 

V / P r 

be solved easily. This was because we put f ^ ~r- z I ( )ctf 

All this was formalism, and its justification depended upon the 
possibility of interpreting the result of the algebraic solution, 
which interpretation could usually, but not always be made. 

There is, however, another method of solution which has the 
same property of reducing the differential equation to an algebraic 

one, in a very similar way. The Laplace transform is defined by: 
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It defines a functional transformation, and relates a time 
function f (t) to a function of the complex variable S which for 
the moment has no physical significance. Examples are: 
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Notice the great similarity - there is one more power of P 



than of S in each case. This is true in general, and is in fact 
the relation by means of which mathematicians have been able to just- 
ify the Heaviside rules for the solution of differential equations. 
The operations of differentiation and integration are analyzed as 



follows : 
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The initial conditions ftBV* are introduced automatically. It 
is seen that the correspondence with the Heaviside method is kept 
here, differentiation corresponding to multiplication of the time 
function by P in the Heaviside system, and to multiplication of 



the complex function by S in the Laplace method. 

The reduction of the differential equation to an algebraic 
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one is then done as follows: 



L ~+ - e,(t) 

Taking the ^transform of both sides, we have: 

+ fci + J-JW] r/M EM (f ) 

i— * 

Or, since the x. transform is linear: 

Ljg'tlM'mM-eW (,t> 

Using our formulas (6) and ( 7 ), we have: 

l [si res) - ?„] +&IC.S) + c -j --Els) 0>) 

This is now reduced to an algebraic ecuation in S, Note that 
the initial conditions give the correct quantities to specify, 
namely the initial current i» an L, and the initial charge on a 
C. Quantities such as the initial voltage on L are not per- 
missible initial conditions, because they can be changed instan- 
taneously. The Heaviside method automatically assumes zero initial 
conditions, because the operational interpretation of p: Jk P t 
requires the (-0) as the lower limit, in order that P and l/P be 
exactly inverse operations, whereas the Laplace transforms have 
(•#G) as the lower limit. 



If we solve ( 11 ) for I (3), we have: 

. e (S) + Li, - Crf 

ICs) - c s 

L S + <K t 

The denominator of ( 12 ) looks very much like the standard 
expression for the steady-state impedance, if we replace S by j (&. 
It can be shown that this is true in general; in other words: 

The relation between the transforms I (3) -and E(S) is 




U) 



u 
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the, steady-state case, if S=^ — - -i-o najjT.rg,' ’ — 

This rule is of great convenience, and enables us to by- 
pass the differential equation . stage unless there is seme 
difficulty with initial, conditions. If initial conditions 
are ignored, the solution assumes them to be zero. 




Third Lecture 



may apply this rule to the series-peaked video amplifier 



circuit discussed in lecture one, as follows: 
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The current applied was a unit step of amplitude 2 , -ass -a-tf IT( 5 ) :l 
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Consulting the^ tables, we find pair 1.109 to be closest, but 
it requires factoring the denominator of the equation: 

rco . k Lli il — 

This requires solving the following equation for S: 



S L S + Lc * 0 



The quadratic formula is used, 
The roots are then: 
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We then can identify: 
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The solution is then: 
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This reduces to: 
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This may be rearranged in the more symmetrical form: 

w =^.[,-^-4 .^ sylA ^ r ' 




which agrees with the result arrived at in the classical method 
in lecture one. 
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Correspondence between Laplace and Fourier transforms may 

Ct/ 0 ° 
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The last formulae are called Fourier transforms. The ex- 



tension of the lower limit to (-©o}Ln the direct transform is made 
easily if t<0 . This is the main difference between 

the Fourier and Laplace transforms, and we will investigate 
further to see just how the complex inversion formula leads to 
time functions vanishing before t~ o 

As a tool useful for later work, let us derive a preliminary 
result. Consider the function 
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crowd together as to increases. However, the area under the curve 



remains^ finite and independent of id, since: 
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Now if this function is multiplied hy a function f (tt) and in- 
tegrated, the rapidly alternating plus and minus values of (sin tot)/t 
will cancel out contributions except in the neighborhood of t equals 
0, where one infinite loop is not cancelled out by another one. 

This means that it is only values of f (t) near to t equals 0 that 
can contribute appreciably to the integral. In the limit, the 
integral is given simply by the value of f (o) multiplied by the 
area under the oscillating curve, which is“i7 : 
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Generalizing slightly by shifting the origin of the oscillating 



curve, we have a formula : 
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Now we return to the Laplace formulae: 
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If we substitute (7) into (8), we have: 
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The order of integration may be reversed, giving 
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Evaluating the inner integral, we have: 
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We then have: 
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This is seen to be the same as (6a) except for the limits of 
integration. From the discussion preceding (6a) it is clear that 
negative values of t must lead to the value f (t) equals 0, since 
contributions to the integral arise only when t and "cross over". 
It is seen that the extra "converging power" of the Laplace 



transform due to the factor & 



_ <r t 



is obtained at the expense of 



cutting off all time functions before t equals 0. If we wish to 
keep F (t) for all values of time, we cannot use the term £ 
but must put S equals J as in equations (2) and (3), which are 
standard Fourier transforms. In general, whatever lower limit was 
taken in the integral defining the direct transformation is the 
time before which the inversion integral gives f (t) equals 0. 



We may also see the above result from the theory of residues. 
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The G^uchy residue theorem is: 
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where the integral encloses the point Z = , and p(z) has no 

pole at Z » <X . The loop integral around any region is then 
equal to £77/ ' ^ residues of all poles enclosed in the region. 
for- t>o , we can define a loop enclosing all poles of the 
function F (S) as below: 




As the integral around the semi- 
circle vanishes for (%s> , the 

integral along the straight line 
is equal to £77/ 

For t < ° . the factor ^ be- 



comes divergent on the semicircle, 
and we have to complete the loop 
by a semicircle to the right of the 
straight line path. Since the new 
loop encloses no poles, we have 



f (t) z o t O . 
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Development of Fourier Integral from Fourier Series 
The usual form of a Fourier series is: 
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This can he put into the complex form by noting that the sine and 
cosine terms of any frequency may be combined as follows: 
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We then have: 
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The similarity between these and the Fourier transform 
formulae is evident. We will not go through the limit process 
because it is too difficult, but it can be seen that a change of 
variable results in the formula 







If now we let o , the limits on variation of t increase without 

limit, and we become interested in much larger values of/), such 
that the product n u^ P remains finite, and equal to to. Then we 
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The other formula becomes: 
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This suggests the definition of an integral, which becomes, 
in the limit of ^ O • 
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so that both of the Fourier integral formulae can be deduced from 
the Fourier series formulae. 

The pfoyiiBiijaal interpretation of a Fourier series is that of ro- 
tating vectors of various frequencies: 






as 0 , this merges into a continuous curve, the projection of 

the end of which is the voltage. 
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In general, Laplace transforms will solve definite problems 
with known initial conditions more easily, while the Fourier trans 
form is best for general arguments about phase and amplitudes, 
since they are directly measurable. 
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Sixth Lecture 

There are certain general conclusions that can he drawn from 
the Fourier analysis that are useful in a number of practical 
cases. For instance, if f(t) and fC*>) are transforms of each 
other, we find that 0 and F ( <-9 e' ; ^ are also trans- 

forms, so that an addition of a linear phase term corresponds to 

S)cfc of tbt. /’haje. , 

a time delay of one second for every cycle / cycle ap u L e u f 
:] lirrrny . Similarly, fit)* and are transforms. 

The time function represents a carrier of frequency fi^ 0 ,mudrriutod ^ 
by f(t) , and the term shows that the spectrum is simply 

that of Ht) but centered around instead of zero. 
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Using these theorems, we can find the effect of suddenly im- 

f re <pne vtc.l<-S 

pressing a sine wave of frequency on a filter which pas ses/f from 
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jdl to ctf2. The impressed is represented as 0<c ] J<lt) e. 

The transform of the bracket is: 
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The network function is given by: 
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The response is then: 
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The integral over negative frequencies may he neglected if 4$ 
is not too far from the positive pass band, so we have: 
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Definitions : 
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The cos. integral may be split into four integrals: 
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If X]_ and x 2 are positive, : r| xj., x 2 r J-X 2 J., and integrals 



(1) and (3) vanish, leaving: 
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If xi < 0 < X2 , we have 
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the 1st and 3rd integrals cancel, since cask? is an odd 
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function. We are left with- 
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the total response is then: 
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if COo is outside the pass band. If 4to is in the band, -we have: 
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Paradox: When the impressed frequency is not in the passband, 

the response ultimately dies out. Then if the impressed frequency 
is suddenly, cut off, the principle of superposition requires that 
another response of reversed sign shall take place. But the net- 
work has physically no driving force at that time, so this response 
must be a frbe transient, with modes characteristic of the network 
a lone . How is it then that Ctb can appear in the solution] 5 j*r 
paatefeawfacgy Sow does the network know what frequency to oBcillate 
at when the drive is removed? How can the network "remember" the 
frequency Wb? 
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Concluding Remarks on Faurier Integrals 

Certain general features of Fourier integrals can often 
be used to get qualitative answers from physical reasoning without 
any equations. Thus, discontinuities in either f(t) or F(w) 
produce wiggles in the other function. This can be seen from the 
"Snake in the Complex plane" interpretation of the Fourier integral. (f>. I f 
If the snake is cut off at a certain frequency, the severed end is 
rotating at that frequency, and wiggles of that frequency appear 
in the time function. This fact may be used in antenna design 
where it is desired to eliminate side lobes. In a broadside array, 
the angular radiation pattern is the Fourier transform of the ampli- 
tude distribution pattern across the array. Since side lobes 
represent wiggles on the 9 pattern, it is seen that the way to 
avoid side lobes is to make the amplitude X distribution fall off 
as smoothly as possible, the smoothest way being a Gauss error 
curve. The 9 pattern is then another Gauss error curve. 

m a? o° 

The way the frequency function falls off as my be 
found by inspection of the function f(t) and the following rules: 



- 1 - 







If F(t) has a: 


F(w) goes down as: 
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No discontinuity in any 
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infinite discontinuity 
{discont. in^fdt) 


F{«») const. 
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infinite doublet disc. 

(discont, injj^dt 2 ) 


F(«f) — 


discont. injj' , *fdt n F(«#) 



rJuL Feason for this is in the varying quality of cancellation in 
the integral defining F(«) , due to the smoothness of f{t). The 
smoother this function, the more perfect the cancellation xkjc&kx 
will be between adjacent half-cycles of the kernel c 
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Analysis of Transmission Lines by the gtfPtft a e Transform 
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Fundamental equations are: 
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taking the Laplace transforms of these equations: 
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in most practical cases we have i = o, vo = 0, so we may set up a 
differential equation for V or I as follows: 
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where , ■-•* 
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The solution of (3) is obviously 
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if we have finite initial values, we may put 
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and we then add to the solution for V the terns: 
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Equations ( 5 ) are recognized as traveling waves, and their 
similarity with the well-know^ steady-state equations suggests 
that we can simply use the steady- state equations in all mani- 
pulations, and then take the inverse transform of the result to 
get our transient solution* 
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Example: Distortion of a Unit Step After Traveling thru a 

Length of Transmission Line. 



Initial voltage = v(o,t) » u(t} 
voltage at point x « v(xjt) 

Transform of initial voltage a VC 0 ) 2 ) ~ T 



transform at point x is 
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attenuation constant 
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so the pair becomes: 



- e 



<T to - ft r* -i 

e 1, l< r vF C7 t?-] 



ViF-tl 



(to< 



the quantity e“^ s + P) to has the inverse transform e «-,(£ ~^e) 
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a "unit impulse" at t : to with area 
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we have: 



This is the same as 
integration with respect to t. 
Therefore, we have: 
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This is a unit step of height e at t = to, plus 

distortion terms due to the integral. The modified Bessel 
Function I| bears somewhat the same relation to an ordinary 
Bessel Function Jn{2) as hyperbolic functions do to trigonometric 
ones. Thus: 
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which corresponds to 
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this analogy can be seen also from the equations 
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Infinite series for In(2) is: 
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the series for Jn(Z) is the same except that signs alternate -f- and - -5 
again corresponding to the relation between hyperbolic and trigo- 
nometric functions. Using the infinite series representation for 



Il(2) , we have 
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The result is expressed as an infinite series.' The accuracy attained 
depends on the number of terms carried. If we carry the first two, 
with m*0 and m» 1, the result will be valid for ^ <C / 



If m : 0 , we have 
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For the voltage 
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we then have, to this approximation: 
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The final value at nay be found from "the following theorem ; 

f(~) -- ■&; * FO) -- e- x 

so that the final distribution is a steady voltage, decreasing 
exponentially with distance from the generator. 




The pulse shape taking only the 1st term of the infinite series is 
given by: 



This predicts a final value of 
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whereas the exact value has been found equal to: 
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This may be approximated by: 
e ' ftb 7 - gr e ~eto[ l ~ rVif-s] 





- f £ 0 o'* U 

c e 9 f» 



This is identical with the value given by the first term of the 
series, showing that the approximate solution is quite good, even for 
large t.±. Incidentally, the equations show that the percentage of 
distortion terms increases as we go down the line, so that eventually 
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the initial sharp rise is lost in distortion terms. The final 
value is always greater than the value of the initial sharp jump. 
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